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Abstract 


The  electronic  structures  of  the  molecules  nitroamlne  and 
nitrosamine  have  been  investigated  by  using  ab  initio  molecular 
orbital  methods.  Double-seta  basis  sets  have  been  used.  Both  planar 
and  non-planar  geoaetries  have  been  considered  and  complete  geoae  try 
optimizations  have  been  performed.  The  results  obtained  are  compared 
with  previous  semi-empirical  and  ab  initio  studies. 
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1.  Introduction 

Nitraaine  and  nitroaaaino  compounds  aro  of  significant  interest  (Feuer 
ot  al.1969)  duo  to  thoir  oxtonsivo  uso  as  essential  ingrodienta  in 
propellante  and  prlaary  and  secondary  explosives.  X-ray  crystallographic 
studies  of  these  compounds  in  conjunction  with  aicrowave  and  electron 
diffraction  studies  indicate  that  these  aolecules  exhibit  a  number  of 
interesting  geoaetrical  features.  For  nitraaine,  apart  froa  the  study  by 
Beevers  et  al.  (1957)  the  aost  exhaustive  study  is  that  of  Tyler  (1963) 
where  the  non-planarity  of  nitraaine  is  shown  by  the  direct  observation  of  the  NH2 
inversion.  However,  there  continue  to  be  different  (Gropen  and  Skancke  1971, 
Ciairaglia  1978)  studies  on  the  contrary,  and  the  fact  that 
diaethyl-nitraaine  is  planar  (Stolevik  and  Radeaacher  1967)  illustrates 
the  necessity  of  further  theoretical  and  experiaental  work  for  the 
accurate  determination  of  the  electronic  structures. 

The  aim  of  the  present  work  is,  therefore,  to  study  the  optimized 
geometries  of  nitraaine  and  nitrosamine  in  both  planar  and  non-planar 
fora,  at  the  Hartree-Fock  level.  The  non-planar  form  of  nitramine  will 
be  taken  from  the  study  of  Tyler  and  the  planar  forms  of  both  the 
aolecules  will  be  taken  from  the  study  of  Gropen  et  al  (1971).  Our 
results  will  be  compared  with  a  number  of  semi-empirical  studies 
(Farainer  and  Webb  1975,  Harris  1973,  White,  Colton,  Lee  and  Rabalais 
1975). 


2.  Detail*  of  Calculation* 


A*  i*  vail  known,  tha  ab  Initio  molecular-orbital  method  predict* 
geometries  with  surprising  accuracy,  in  spite  of  the  neglect  of 
correlation  energy.  However,  there  appears  to  be  only  two  ab  initio 
studies  in  the  literature  for  these  molecules.  The  study  by  Gropen 
et  al.  (1971)  considered  only  planar  geometry  and  no  geometry 
optimization  was  performed.  The  study  by  Duke  (1978)  used  minimum 
basis  sets  and  only  partial  geometry  optimization  was  performed.  The 
NC>2  and  NH2  geometries  were  kept  fixed;  only  the  H-H  distance  and  the 
out-of-plane  angle  of  the  NH2  group  were  varied. 

To  improve  these  results,  we  have  used  a  double-zeta  basis  set, 
namely  the  3-21G  basis  set  (Binkley,  Pople  and  Hehre  1980).  This  basis 
consists  of  an  s-type  inner  shell  function  with  3  Gaussians,  an  inner 
set  of  valence  a-  and  p-  type  functions  with  2  Gaussians  and  another 
outer  sp  set  with  1  Gaussian.  This  basis  should  give  reasonable  bond 
angles  and  bond  lengths,  which  are  known  to  be  overestimated  in  a 
minimal  basis-set  study.  Geometry  optimization  was  performed  by  using 
the  method  of  Schlegel  (1982).  This  optimization  algorithm  evaluates 
and  utilizes  the  gradients  each  time  the  energy  is  computed  and  the 
second  derivative  matrix  is  updated.  This  is  followed  by  a  one¬ 
dimensional  search  using  the  second  derivative  matrix. 


3.  Results  and  Discussion 

The  results  obtained  axe  presented  in  Tables  1  -  12. 

In  Tables  1  and  2,  the  non-planar  experimental  geometries  of 
nitraaine  and  the  optimized  geoaetry  at  the  Hartree-Fock  3-2 1G  level 
are  presented.  The  N-N  distance  obtained,  1.472  A,  is,  not  sur¬ 
prisingly,  too  large.  This  is  to  be  compared  with  the  N-N  distance, 
1.354  A,  obtained  at  planar  geoaetry.  Obviously,  there  is  a  strong 
tendency  to  planarity  at  shorter  N-N  distances.  MINDO  predicts 
nitraaine  to  be  a  planar  molecule,  MIND02  gives  1.21  A  for  N-N 
distance,  while  MINDO 3  gives  1.29  A. 

In  tables  3  and  5,  we  present  the  calculated  orbital  energies  for 
nitraaine  in  non-planar  and  planar  forms.  The  nature  of  each  orbital 
is  indicated  in  terns  of  w-bonds  and  aore  or  less  localized  o  bonds. 

The  four  lowest  occupied  orbitals  are  the  inner  shell  s-orbitals  on  the 
four  heavy  atoas.  This  feature  is  also  supported  by  the  MIND02  and 
MIND03  orbital  energies.  The  photoelectron  spectrua  of  nitraaine 
Indicates  that  the  two  highest  occupied  levels  are  of  n  symmetry.  Our 
calculations  indicate  that  while  this  is  true  for  planar  nitraaine,  in 
non-planar  nitraaine,  only  the  highest  occupied  level  is  of  «  symmetry. 
Also,  in  the  study  by  Gropen  et  al.,  one  lone  pair  is  found  between  the 
two  o  (N-H)  orbitals  and  the  lowest  r  orbital  was  found  to  be  almost 
degenerate  with  the  highest  of  these.  We  do  not  observe  such 
degeneracy  and  thus  it  seems  unlikely  that  the  lowest  w  orbital  is 
alaost  completely  delocalized  as  has  been  claiaed  in  the  literature. 
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In  Tables  6  and  7,  the  optimized  planar  geometry  and  the  orbital 
energies  for  nitrosaalne  are  presented.  Here  the  three  lowest  orbitals 
are  the  inner  shell  s-orbitals  on  the  three  heavy  atoms.  A  study  of 
population  analysis  indicate  coaplete  agreeaent  with  the  conclusions 
reached  by  Gropen  et  al.:  the  three  highest  occupied  9  orbitals  are 
responsible  for  the  lone  pairs  on  nitrogen  and  oxygen.  Also,  the  two 
lowest,  5a'  and  4a',  are  aainly  responsible  for  the  9  (N-0)  and  9  (N-N) 
bonds . 

Table  8  contains  the  results  of  the  total  energies,  before  and 
after  optimization.  Since  nitraaine  in  planar  fora  has  clearly  lower 
energy  than  nitraaine  in  non-planar  fora  and  since  experiaental  study 
by  Tyler  (1963)  indicate  results  on  the  contrary,  further  experiaental 
research  in  this  direction  could  be  very  fruitful. 

Table  9,  10,  and  11  present  the  results  of  gross  atomic  population 
study.  In  the  non-planar  fora,  the  nltro  nitrogen  atom  is  positively 
charged,  naaely  +0.51  while  the  amine  nitrogen  atom  is  negatively 
charged,  naaely  -0.51.  This  is  from  Mulllken  population  analysis. 
Results  from  Lowdin  population  analysis  indicate  that  the  nitro 
nitrogen  atoa  is  positively  charged,  +0.46  but  the  amine  nitrogen  is 
strongly  negatively  charged,  namely  -0.71.  The  same  qualitative  trends 
are  predicted  from  the  CND0  charges  which  are  in  reasonable  agreement 
with  the  IND0  charges  reported  by  White  (1975).  It  is  interesting  to 
note  that  for  nitramlne  in  planar  form,  the  conclusions  are  completely 
opposite.  We  note  as  a  sideline  that  MINDO  methods  do  not  work  very 
well  in  the  study  of  the  population  charges  of  these  molecules, 
probably  because  these  molecules  are  very  different  from  those  used  in 


the  paraaeterisation  of  ouch  methods.  It  is  also  to  be  remembered  that 
the  predicted  populations  in  any  ab  initio  calculation  are  very 
sensitive  to  the  choice  of  basis  set. 

Finally  in  Table  12,  we  compare  our  results  for  the  dipole  moment 
of  nitraaine  with  the  available  results  in  the  literature.  In  the 
non-planar  fora,  our  value  for  the  dipole  nonent  is  the  closest  to  the 
experimental  value.  However,  since  total  energy  is  considered  aore 
reliable  in  ab  initio  aethod,  nitraaine  in  planar  fora  still  appears  to 
be  aore  stable  theoretically.  Again,  further  experimental  work  would 
be  very  welcome. 
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TABLE  1 

Experimental  bond  lengths  and  bond  angles  for  nitramine  (lengths  in.nn) 


Beevers  and 

Trotman-Dickinson  (1957) 

X-ray 

Tyler  (1963) 

Microwave 

N-H 

0.1005  i  0.001 

N-H 

0.140 

0.1427  ±  0.0002 

N-0 

0.118 

0.1206  (assumed) 

0N0 

130*8'  ♦  15’ 

HNH 

115*11’  i  2* 

out  of  plane 

51*47’  ♦  1* 

TABLE  2 

Optiaized  Non-planar  Geoaetry  of  Nitraaine  Coordinates  are  in  Angst^oaa. 


Atoa 

X 

y 

z 

H(  1 ) 

0.0145 

-0.0021 

0.0 

N(2) 

1.4845 

0.0754 

0.0 

0(1) 

-0.5963 

1.0803 

0.0 

0(2) 

-0.4676 

-1.1209 

0.0 

H(l) 

1.7172 

0.6055 

0.8295 

H(2) 

1.7172 

0.6055 

-0.8295 
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TABLE  3 

Orbital  energies  for  nitroaaine  in  a.  u. 


Orbital 

Orbital  energy 

Nature  of  Orbital 

1  (la*) 

-20.531190 

Is 

2  (2a‘) 

-20.511603 

Is 

3  (3a‘) 

-15.810688 

Is 

4  <4a') 

-15.574496 

Is 

5  (5a*) 

-  1.631171 

0 

6  (6a') 

-  1.424051 

0 

7  (7a1) 

-  1.261226 

0 

8  (8  a’) 

-  0.905639 

0 

9  (la") 

-  0.786352 

» 

10  (9a') 

-  0.775029 

0 

11  (10a') 

-  0.723289 

0 

12  (2a") 

-  0.665820 

* 

13  (11a') 

-  0.538030 

a 

14  (12a’) 

-  0.499191 

a 

15  (13a1) 

-  0.467081 

a 

16  (3a"> 

-  0.443007 

i 

TABLE  4 


Initial  and  Optimized  Planar  Geoaatri.es  of  Nitraaine  Coordinates  are  in 
Angetroaa. 


Atom 

xin 

*in 

X  a, 

opt 

yopt 

H(l) 

0.0 

0.0 

0.0 

-0.0065 

N(2) 

0.0 

1.3816 

0.0 

1.3476 

0(1) 

1.0926 

1.8911 

1.1170 

1.9019 

0(2) 

-1.0926 

1.8911 

-1.1170 

1.9019 

H(l) 

0.8474 

-0.5398 

0.8937 

-0.4426 

H(2) 

-0.8474 

-0.5398 

-0.8937 

-0.4426 
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TABLE  3 


Optimized  Orbital  Energies  for  nitramlne 

in  planar  form  in  a.u. 

Orbital 

Orbital  energy 

Nature  of  Orbital 

1  (la*) 

-20.495737 

la 

2  (2a*) 

-20.495366 

la 

3  (3a') 

-15.814295 

la 

4  (4a') 

-15.580724 

la 

5  (5a’) 

-  1.617184 

e 

6  (6a’) 

-  1.401714 

9 

7  (7a') 

-  1.282449 

9 

00 

00 

-  0.924754 

9 

9  (9a') 

-  0.799121 

9 

10  (10a’) 

-  0.769095 

9 

11  (la") 

-  0.746520 

V 

12  (11a’) 

-  0.671962 

9 

13  (12a') 

-  0.497412 

9 

14  (13a') 

-  0.482884 

9 

15  (2a") 

-  0.472506 

9 

16  (3a") 

-  0.437638 

9 

u 


TABLE  6 

Initial  and  Optiaized  Planar  Geoaetriee  of  Nitroeaaine.  Coordinate*  are  in 
Angatroae. 


Atoa 

xin 

yin 

xopt 

yopt 

N(l) 

0.0 

0.0 

0.0060 

-0.0062 

M(2) 

0.0 

1.3436 

-0.0164 

1.3413 

0(1) 

1.0926 

1.8531 

1.1001 

1.8513 

HU) 

0.8474 

-0.5398 

0.8896 

-0.4737 

H(2) 

-0.8474 

-0.5398 

-0.8690 

-0.4746 
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TABLE  7 

Optiaized  Orbital  Energies  for  nitrosaaine  in  a.  u. 


Orbital 

Orbital  Energy 

Nature  of  Orbital 

1  (la') 

-20.489050 

Is 

2  (2a') 

-15.659971 

Is 

3  (3a*) 

-15.546273 

Is 

4  (4a') 

-  1.523223 

0 

5  (5a') 

-  1.248262 

0 

6  <6a') 

-  0.915149 

0 

7  (7a') 

-  0.770635 

a 

8  (8a') 

-  0.696589 

9 

9  (la") 

-  0.623923 

9 

10  (9a’) 

-  0.586597 

a 

11  (2a") 

-  0.424889 

* 

-  0.414389 


a 


TABLE  8 

Total  oloctronlc  energies.  Values  In  a.  u. 


Nolacula  Geoaetry  Total  Enorgy 


Mltraalne 


Mltraalne 


Nitrosaalno 


Mon-Planar  (Expt.) 
Mon- Planar  (Opt.) 
Planar  (Thao.) 
Planar  (Opt.) 
Planar  (Thao.) 
Planar  (Opt.) 


-258.004191 

-258.118307 

-258.128685 

-258.137938 

-183.764893 


-183.768075 


TABLE  9 

Groat  atoaic  population  for  nitraaina  in  non-planar  geoaatry. 


Total  Charga 
(Lovdin) 


Atoa 


Total  Charga 
(Mullikan) 


■(1)  6.49297  6.54222 
11(2)  7.51055  7.29311 
0(1)  8.39306  8.33289 
0(2)  8.33329  8.24663 
H(l)  0.63507  0.79258 
H(2) 


0.63507 


0.79258 


TABU  10 

Gro»a  atoaic  population  for  nltraaina  in  planar  geoaatry. 


Atoa  Total  Charga 

(Nulllkan) 


N(l) 

N(2) 

0(1) 

0(2) 

H(l) 

H(2) 


7.55396 

6.45242 

8.40378 

8.40378 

0.59033 


Total  Charge 
(Lovdln) 


7.29915 

6.48963 

8.34106 

8.34106 

0.76455 


0.59033 


0.76455 


TABLE  11 

Grots  atomic  population  for  nitrosaaine  in  planar  geometry. 


Atom 


Total  Charge  Total  Charge 

(Mulliken)  (Lowdin) 


N(l)  7.62222  7.31078 
N(2)  6.71360  6.81728 
0(1)  8.41915  8.29476 
H(l)  0.62775  0.79472 
H(2) 


0.61729 


0.78245 


TABLE  12 


Coaparleon  of  dipole  aoaents  of  nitraaine. 


Method 

Nitraaine 

CNDO  (Duke  1978) 

4.44 

MIND03  (Duke  1978) 

3.95 

INDO  (Whitman  and  Hornback 

4.12 

1969) 

Ab  Initio  (Gropen  et  al.  1971) 

4.42 

Ab.  initio  (Duke  1978) 

3.75 

This  work 

(non-planar) 

3.62 

This  work 

(planar) 

4.66 
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ABSTRACT 


The  recent  availability  of  artificial  intelligence 
programs  and  Lisp  machines  brings  a  completely  new  class 
of  computing  tools  to  the  physicist.  In  this  paper  we 
try  to  provide  a  pedagogical  introduction  to  this  not  -o 
common  class  of  computing  by  demonstrating  some  of  the 
general  applications  of  these  programs,  by  showing  one 
of  its  applications  in  solid  state  physics  research,  and 
by  suggesting  other  research  and  classroom  applications. 


I.  INTRODUCTION 


During  the  past  decade  the  use  of  digital  computers  has  be 
come  an  Integral  part  of  research  in  solid  state  physics.  In 
most  cases  the  uses  of  computers  have  been  primarily  either  for 
numerical  calculations  or  for  control,  data  acquisition,  and 
data  analysis  in  experiments.  In  this  paper  we  try  to  point 
out  a  different  type  of  use,  namely  for  algebraic  (or  symbolic) 
manipulation.  This  type  of  computer  programming  is  part  of  a 
larger  field  which  Is  known  as  artificial  intelligence  and  Is 
expected  to  make  a  very  large  impact  in  the  next  generation  of 
computers. 

Although  a  munber  of  researchers  in  physics  have  used 

computers  for  algebraic  programming,  it  is  still  an  unfamiliar 

topic  to  a  large  fraction  of  physics  faculty  and  students. 

Thus  In  the  first  part  of  this  paper  we  take  a  pedagogical 

approach  and  introduce  the  use  of  algebraic  programming  via 

examples.  In  subsequent  parts  we  present  its  application  to  a 

particular  problem  in  solid  state  physics,  namely  the  formu- 

lation  of  band  structure  calculations  based  upon  the  linear 

combination  of  atomic  orbitals  (LCAO)  method.  The  examples 

presented  in  the  rest  of  the  paper  are  based  on  the  "computer 

language"  REDUCE.1  The  easy  availability  of  this  language 

2 

over  the  more  powerful  ones,  like  MACSYMA,  was  the  reason  for 
this  choice. 


II.  AL6EBRAIC  PROGRAMMING 


The  primary  purpose  of  algebraic  manipulation  programs  Is 
to  process  formal  mathematical  expressions,  without  any  particular 
concern  for  their  numerical  values.  Basically  these  expressions 
are  processed  In  a  fashion  similar  to  the  one  used  In  an  algebra 
or  calculus  class.  Since  almost  all  of  the  readers  of  this  paper 
are  familiar  with  numerical  programming  using  computer  languages 
like  BASIC  or  FORTRAN,  we  think  It  Is  most  Instructive  to  present 
an  example  comparing  the  algebraic- programming  method  with  the 
numerical  programming  method. 

Example  I 

Let  us  consider  the  expression 

(X2  -  Y2)  /  (X+Y).  (1) 

If  we  write  a  FORTRAN  program  as 

Z  =  (X**2  -  Y**2)  /  (X+Y)  (2) 

PRINT  (5,*)  Z  (3) 

END  (4) 

and  execute  it.  the  result  will  be  completely  meaningless  and 
some  compilers  will  warn  that  the  X  and  Y  variables  are  "un> 
defined."  Running  a  similar  program  from  a  BASIC  interpreter  will 
produce  a  similar  error  message.  If  we  run  a  program  in  REDUCE  as 

Z  :  = 


(X**2  -  Y**2)  /  (X+Y); 


(5) 


the  result  will  be  (X  -  Y). 

The  Min  difference  between  the  two  types  of  programing  is 
Illustrated  In  tbe  above  example.  In  the  numerical  programing 
method,  one  must  assign  a  numerical  value  to  each  variable  before 
an  expression  containing  It  is  evaluated,  whereas  in  algebraic 
programing  one  obtains  an  answer  which  Is  true  for  all  different 
values  of  the  variables.  Since  in  solid  state  physics  problems, 
as  well  as  in  a  number  of  other  branches  of  science,  one  needs  to 
translate  a  physical  model  Into  mathematical  expressions  before 
starting  numerical  programing,  in  most  of  these  cases  the 
algebraic  programing  method  could  be  used  to  obtain  those  ex¬ 
pressions.  In  table  I  we  have  presented  a  few  examples  using 
REDUCE.  On  one  line  the  operation  is  defined,  followed  by  a  semi¬ 
colon.  On  the  lines  after  the  semicolon  the  results  from  the 
execution  of  REDUCE  are  printed.  Examples  of  differentiation 
and  integration  are  -given.  The  statements  starting  with  ON  or 
OFF  are  control  directives  to  the  REDUCE  program  and  may  be 
ignored  during  the  reading  of  this  paper  without  loss  of  under¬ 
standing.  The  computer  used  for  these  examples  was  a  DEC  2060 
running  with  the  T0PS-20  operating  system,  but  these  examples 
can  be  run  in  any  implementation  of  REDUCE  yielding  the  same 
results. 


III.  LCAO  BANDS TRUCTURE 


In  this  section,  as  a  specific  application,  we  present  the 
use  of  REDUCE  In  obtaining  the  expressions  for  the  linear  com- 
blnatlon  of  atomic  orbitals  (.LCAOl  Hamiltonian  matrix  elements 
for  electrons  In  a  crystal.  An  Introduction  to  the  band  theory 
In  general  and  the  LCAO  method  In  particular  can  be  obtained 

from  many  solid  state  physics  texts.  Also,  the  original  paper 
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by  Slater  and  Koster.  and  a  review  article  by  Nussbaum.  con¬ 
tain  a  comprehensive  dtscusslon  of  the  theory  Involved,  in 
light  of  this,  a  detailed  Introduction  will  be  omitted  here. 

This  method  has  been  a  very  effective  starting  point  for  the 

6  7  8 

study  of  defects,  surfaces,  the  electron-phonon  Interaction, 

Q 

and  Fermi-surface  properties. 

In  the  LCAO  formulation  of  the  band  structure  problem,  one 
expresses  the  conduction  electron  wave  function  as  a 

linear  combination  of  atomic  orbitals  (hence  the  name  LCAO): 


—  z 
/N  mj 


ei*,Rmd 


(6) 


where  ♦jm(r)  Is  the  atomic  orbital  centered  at  site  m  and  the 
summation  j  runs  over  all  the  atomic  orbitals  considered,  e.g. 
orbitals  of  symmetry  s,  x,  y,  xy,  yz,  etc.  The  energy  levels 
and  wave  functions  for  the  electrons  are  obtained  by  diagonal¬ 
izing  the  matrix  H^,  where 
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» 


H,J (7) 

and  H  ts  the  Hamiltonian,  The  above  expression  for  after 
a  commonly-used  approximation  called  the  two  center  approxi¬ 
mation4,5  can  be  written  as 

p (.neighbors *PEij(V  (8) 

where  ftp  represents  the  lattice  sites  and  EtJ  Is  a  unique 
function  of  1,  j,  and  ifp,  The  functional  dependence  of  Ejj  on 
ltp  can  be  further  subdivided  Into  radial  and  angular  parts. 

The  E|j  expressions  for  a  few  different  values  of  1  and  j  are: 

Esx  "  *(sP°>  (9) 

E  2  2  -  itf/2U2-m2}(sdo}  (10) 

s.x  -y 

E  2  7  "  «-(n2-( A2+m2)/2)CpdCT)  -  /Imn2(pdir)  (11) 
x,3z  -r 

where  l,  m,  and  n  are  the  direction  cosines  of  7tpi  representing 
the  angular  part,  and  (sso),  ( pp-rr) ,  etc.  are  two-centered  bond 
parameters,  representing  the  radial  part.  A  complete  list  of 
the  functions  is  given  in  references  4  and  5.  A  general 
procedure,  Which  can  also  be  implemented  using  REDUCE,  is 
given  by  Sharma.10 

In  evaluating  H^,  the  contribution  from  the  radial  part 
Is  usually  determined  from  numerical  calculations  or  from  a 


i 
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parametric  fit,  and  that  from  the  angular  part  Is  obtained  by 
analytically  sunlng  over  the  ftp.  As  an  example,  we  will  con- 
stder  the  contribution  to  for  the  nearest  neighbors  In  a 
body-centered  cubic  (.bcc)  lattice.  The  number  of  nearest 
neighbors  In  this  case  Is  eight  and  their  coordinates  and 
(t  m  n)  values  are  listed  In  Table  II.  Substituting  these 
values  of  U  *  n)  and  the  form  of  from  equation  9  into 
equation  8,  one  gets 

-  21  s1n(kxa)(spo)  .  (12) 

where  a  Is  the  lattice  constant.  This  expression  Is  the 
contribution  of  the  nearest  neighbors  In  a  bcc  lattice  to  H^. 
Although  the  contributions  from  other  matrix  elements,  and  for 
different  lattices  can  be  obtained  in  an  identical  fashion,  the 
process  Is  extremely  tedious,  error-prone,  and  time  consuming. 
Table  III  contains  a  program  in  REDUCE  which  evaluates  the  sum 
over  ftp,  given  the  ftp  vectors.  Readers  interested  in  getting  a 
feel  for  the  advantage  of  using  this  program  may  try  evaluating 
one  of  the  d-d  matrix  elements,  or  may  look  in  Ref,  4  or  5  at 
the  complicated  results  for  s,  p,  and  d  orbital  symmetries. 
Readers  not  .familiar  with  the  programming  language  REDUCE  should 
only  be  concerned  with  the  remarkable  brevity  of  a  program 
which  yields  the  complicated  LCAO  expressions  for  arbitrary 
angular  momentum  orbital  basis  functions. 
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IV.  OTHER  APPLICATIONS 

Artificial  In talll genet  programing  has  countless  appll- 
cations.  These  applications  Include  compl ica ted  Integrations, 
differentiations.  Halt  evaluations,  suaaatlons,  and  application 
of  recursion  relations  to  give  analytic  expressions.  These 
prograas  are  very  efficient  In  performing  the  more  tedious 
kinds  of  operations  In  which. hunans  have  a  high  rate  of  error. 

Since  most  students  tn  physics  as  well  as  the  other  sciences 
are  acquainted  with  personal  computers  and  mainframe  computers, 
and  use  them  In  both  graduate  and  undergraduate  studies,  an 
Introduction  to  artificial  Intelligence  programming  at  an  early 
stage  In  the  student's  career  is  now  appropriate.  This  can  be 
done  In  many  physics  classes,  but  certainly  in  courses  in 
mathematical  methods  of  physics  at  any  level. 

A  few  other  examples  of  useful  REDUCE  programs  are  programs 
which  produce  explicit  expressions  for  the  following: 

1)  orthogonal  polynomials  of  arbitrary  order 

2)  spherical  harmonics  of  arbitrary  1  and  m 

3)  spherical  Bessel  functions  of  arbitrary  order 

4)  Clebsch-Gordon  coefficients. 

A  comment  on  the  last  item  is  informative  and  illustrative  of  the 
power  of  REDUCE.  Many  raan-years  of  effort  have  been  expended  in 
computing  Clebsch-Gordon  coefficients  of  high  order.  Numerical 
calculations  are  difficult  because  of  factorial  functions  of 
large  arguments.  Extensive  tables  of  the  coefficients  have  been 


published,  both  fn  prise  number  representation  and  decimal 
representation.  Using  REDUCE  we  have  written  a  simple  program 
wbtcb  very  quickly  gives  an  analytic  expression  for  or  a  nu¬ 
merical  value  of  the  coefficient  for  arbitrary  angular  momenta 
without  any  of  the  computational  difficulties  encountered  In 
previous  methods.11 


V .  CONCLUSION 


ife  have  demonstrated  the  utility  of  algebraic  programming 
by  presenting  a  specific  example  of  how  It  has  been  used  In 
solid  state  physics  research. ^  The  advantages  of  this  method 
are  threefold.  First,  the  tedium  of  deriving  many  expressions 
using  the  same  method  is  greatly  reduced.  Second,  a  great  deal 
of  time  can  be  saved  by  ustng  the  computer  to  do  repetitive 
derivations.  Finally,  the  chances  of  error  In  the  derivations 
and  In  transcription  of  the  equations  are  minimized.  The  latter 
Is  due  to  the  fact  that  REDUCE  has  the  capability  of  writing 
the  results  In  FORTRAN, 

With  the  recent  stress  on  arltlflcal  intelligence  and  fifth 
generation  computers,  there  is  a  growing  need  for  more  re¬ 
searchers  and  students  to  become  familiar  with  programming 
languages  like  REDUCE.  Not  only  are  they  useful,  they  may  soon 
become  a  necessary  part  of  any  scientific  researcher's  tools. 
Artificial  intelligence  may  remove  the  bottleneck  between  the 
conception  of  a  solution  to  a  problem  and  its  execution,  giving 
the  student  and  researcher  more  time  to  think  and  learn,  as  well 
as  the  ability  to  attack  problems  of  greater  complexity. 

We  are  convinced  that  in  the  future  the  use  of  artificial 
Intelligence  programs  will  be  just  as  important  a  part  of  training 
for  scientists  as  learning  how  to  read  a  voltmeter  or  use  a 
microscope. 
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Table  Captions 

Table  I.  Examples  of  the  use  of  REDUCE  for  differentiation 
and  integration  of  selected  functions. 

Table  II.  Nearest-neighbor  positions  and  direction  cosines  for 
the  BCC  lattice.  The  nearest-neighbor  distances  have 
been  multiplied  by  2/a,  where  a  Is  the  lattice  con¬ 
stant,  and  the  direction  cosines  have  been  multiplied 

by  47. 

Table  III.  A  REDUCE  program  which  was  used  for  the  evaluation 
of  Slater-Koster  matrix  elements. 


Table  I 


Tht  following  art  examples  of  differentiation: 


..  v- 


;v  .  •••-.■ 


i  •  I  . 


Th*  syntax  ii  DP{  F  ,  x  ,  n), 

Vtfiara  function  F  is  differentiated  n  times  with  respect 
to  the  variable  x  ; 


DP(  X**3  ,  x  ,  1); 


3*X 


DP(  X**H  *  SIN(X)  ,  X  ,  2); 


N 


(X  *(2*COS(X)*X*H  -  SIN(X) *X  ♦  SIN(X)*N  -  SIN(X) *H) ) /Y* 


DP(  X**(X**X)  ,  X  ,  3) ; 


(X  ♦  X)  (2*X)  6  3  (2*X)  5  3  (2*X) 

-  <X  *LOG(X)  *X  ♦  3*x  *LOG(X)  *X  ♦  3*X  * 


4  3  (2*X)  4  2  (2*X)  3  3 

LOG(X)  *X  ♦  3*X  *LOG(X)  *X  ♦  X  *LOG(X)  *X  ♦ 


(2*X)  3  2  ( 2*X)  2  2  (2*X) 

6*X  *LOG(X)  *X  ♦  3*X  *LOG(X)  *X  ♦  3*X  * 


2  (2*X)  (2*X)  X  5 

LOG(X)  *X  +  3*X  *LOG ( X ) *X  ♦  X  ♦  3*X  *LOG(X)  * 


3  X  4  3  X  3  3  X  3 

X  +  9*X  *  LOG  ( X )  *X  *  9*X  *  LOG  ( X )  *X  ♦  12*X  *LOG(X)  * 


2  x  2  3  X  2  2  X  2 

X  +  3*X  *LOG(X)  *X  +•  2 1  *X  *LOG(X)  *X  -  3*X  *LOG(X)  * 


X  2  X  XX 

X  +  9*X  *LOG ( X ) *X  ♦  6*X  *LOG ( X )  *X  ♦  6*X  *X  -  3*X  ♦ 


43  33  23  22 

LOG  ( X )  *X  +  3  *  LOG  ( X  )  *X  ♦  3*LOG(X)  *X  ♦  6*LOG(X)  *X 


3  2  2  3 

♦  LOG(X) *X  +  9* LOG ( X )  *X  -  4*LOG(X)*X  ♦  3*X  +  2 ) ) /X 


DP<  A*X  *  ASIN(X)  ,  X  ,  2); 


2.2  42 

( SQRT ( X  -  1)*A*I*(  -  X  ♦  2 ) ) / ( X  -  2*X  +  1) 


Table  II 

8CC  Lattice 


Nearest-neighbor  position  (x  2/a) 


Direction  cosines  (xJT) 


Comment  Thle  program  derives  the  equation  for  TB  bandstructure; 
OFF  NERO;  OFF  ECHO; 

Comment  Read  in  the  E(l,m,n)s  and  the  R  vectors; 
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Abstract 


A  detailed  study  of  Li^i  Lij  and  Li^  clusters  is  reported  at  the 
ab  initio  level.  Fourteen  different  geometries  of  the  clusters  have  been 
considered.  The  most  stable  geoaetrical  forms  at  the  Hartree-Fock- 
Configuration-Interaction  level  are  found  and  compared  with  previously 
published  theoretical  and  experimental  results. 
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I .  Introduction 

In  recent  years,  there  has  been  widespread  interest  in  theoretical 
investigation  of  aetal  clusters  (Schaefer  1975).  The  understanding  of  the 
properties  of  saall  clusters  is  of  signal  iaportance  in  aany  areas, 
particularly  dispersed  aetal  catalysis.  Lithiua,  as  the  simplest  aetal, 
has  been  the  subject  of  aany  of  these  studies  (Stoll  and  Preuss  1972, 
Marshall,  Blint  and  Kuna  1976,  Borisov  1976,  Gerber  and  Schumacher  1978, 
Heraann  and  Bagus  1978,  Beckmann  at  al.  1979)  primarily  because  siaple 
quantum  cheaical  methods  as  well  as  sophisticated  self-consistent-field 
aethods  can  be  eaployed  for  evident  reasons.  As  a  result,  a  number  of 
results  of  different  quality  on  different  clusters  is  now  available.  For 
exaaple,  ab  initio  Self-Consistent-Field  (Davies  and  Del  Conde  1976), 
Configuration  Interaction  (Kress,  Carberry  and  Kuczynski  1978),  Coupled- 
Electron-Pair-Approximation  (Meyer  1977)  and  Diatomics-In-Molecules 
(Companion,  Steible  and  Starshak  1968)  studies  of  Li  clusters  have  been 
published  but  the  obtained  results  are  hardly  in  agreement  with  each  other. 
In  this  study  and  the  studies  to  follow,  we  want  to  use  one  method  only, 
namely  Hartree-Fock  followed  by  full  configuration  interaction  in  the 
entire  spectrum  of  diatomic  molecule  to  very  large  clusters.  The  ultimate 
aia  is  to  determine  precisely  what  size  and  shape  of  the  cluster  which  can 
be  used,  to  a  reasonable  approximation,  to  simulate  the  effects  of  the 
infinite  solid.  This  work  concentrating  on  small  L^.  Li^  and  Li^  clusters 
is  the  first  one  in  a  series  of  exhaustive  study  of  metal  clusters  at  the 
Hartree-Fock-Configuration-Interaction  level.  It  is  hoped  that  auch  a 
consistent  study  will  help  resolve  the  controversies  surrounding  the  size, 
shape  and  stability  of  saall  Li  clusters. 
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XI.  Geometry  of  Cluster* 

Clusters  investigated  in  this  work  are  taken  froa  different  studies 
available  in  the  literature  and  are  shown  in  Fig.  1.  Triangular, 
equidistant  and  general  Li^  clusters  are  investigated.  The  paraaeters  are 
taken  froa  a  coaparable  CHDO/BW  (Coaplete  Neglect  of  Differential  Overlap) 
study  (Skala  1981)  since  it  is  known  that  such  paraaeterisation  is  able  to 
give  the  quantities  of  priaary  interest,  i.e.  bond  lengths,  bond  angles  and 
binding  energies  to  a  fairly  good  accuracy  (Boyd  and  Whitehead  1981).  Of 
course,  it  is  expected  that  a  HF  followed  by  Cl  calculation  will  improve  the 
results  of  seaieapirical  studies  considerably.  Ten  different  structures  of 
Li^  are  taken  into  account:  tetrahedral  (tet),  triangular  (tri),  square 
(sq),  linear  equidistant  (eq),  linear  consisting  of  two  diatomic  molecules 
(aol),  general  linear  (lin),  oblong  formed  by  two  parallel  diatomic  molecules 
(obi),  T-shaped  (T),  rhombic  formed  by  two  diatomic  molecules  (rho)  and 
parallelopiped  (para). 

HI.  Computational  Method 

One  of  the  primary  considerations  involved  is  determination  of  the 
type  of  basis  set  to  be  used.  Gaussian-type  basis  sets  used  in  ab  initio 
molecular  orbital  computations  usually  involve  some  compromise  between 
computational  cost  and  accuracy.  A  considerable  increase  in  computational 
efficiency  can  be  achieved  if  the  exponents  of  the  Gaussian  primitives  are 
shared  between  different  basis  functions.  At  the  split-valence  level,  this 
is  usually  exploited  by  sharing  primitive  exponents  between  s  and  p 
functions  for  the  valence  functions.  Accordingly,  for  lithium,  we  have 


4 


used  here  the  so-called  3-21  G  basis  set  (Binkley,  Pople  and  Hehre  1980). 
This  basis  consists  of  s-type  inner  shell  function  with  3  Gaussians,  an 
inner  set  of  valence  a-  and  p-type  functions  with  2  Gaussians  and  another 
outer  sp  set  with  1  Gaussian. 

Once  the  basis  set  has  been  decided  upon,  Hartree-Fock  calculations 
have  been  carried  out  followed  by  full  Cl  involving  all  single  and  double 
excitations.  The  use  of  the  method  of  configuration  interaction  is  crucial 
since  it  is  well  known  (Das  1967)  that  conventional  Hartree-Fock  aethods  do 
not  describe  Lin  systeas  accurately.  These  Cl  calculations  have  been 
perforaed  within  the  fraaework  of  the  Graphical  Unitary  Group  Approach 
(GUGA-CI)  of  Brooks  and  Schaefer  (1979).  It  is  to  be  noted  that  GUGA-CI 
approach  is  computationally  well  tractable  compared  to  conventional  Cl 
aethods. 

IV.  Results 

The  results  of  our  calculations  are  presented  in  Tables  1  and  2.  The 
parameters  are  listed  in  Column  2  of  Table  1 .  AE  indicates  the  energy 
separation  between  the  lowest  unoccupied  molecular  orbital  and  the  highest 
occupied  molecular  orbital  and  EtQt  is  the  total  energy,  i.e.  the  core 
repulsion  energy  plus  the  valence  electron  energy.  is  the  binding 
energy  per  atom  given  by 

h  ■  <NILi  *  ELiN),N  (I> 

where  N  is  the  number  of  atoms  in  the  cluster.  Table  2  lists  the  first  and 
the  second  ionization  potentials.  The  first  ionization  potential  is 
obtained  from  Koopman's  theorem 

homo 


I.  P.  -  -E, 


(2) 


and  the  second  ionisation  potential  was  obtained  from  the  orbital  below  the 
highest  occupied  aolecular  orbital. 

For  the  separation  distance  was  taken  to  be  the  equilibrium 
experiaental  separation.  Our  value  for  the  first  ionization  potential  4.89 
eV  compares  favorably  with  the  experiaental  ionisation  potential  4.96  eV. 
The  difference  could  be  attributed  to  the  basis  set.  The  Cl  was  fairly 
snail,  involving  only  553  configurations,  counting  all  single  and  double 
excitations. 

For  Li^,  the  conclusions  agree  with  the  results  of  Skala  (1981)  and 
Kendrick  and  Hillier  (1977)  at  the  Hartree-Fock  level.  This  ioplies  that 
at  this  ab  initio  level,  the  nost  stable  Li^  cluster  has  the  fora  of  a 
non-equilateral  triangle.  The  total  energy  of  linear  equidistant  Li^  is 
found  to  be  lower  than  that  of  the  triangular  Li^.  This  agrees  with  the 
results  of  Skala  (1981)  and  Coapanion  (1978)  but  in  contrast  to  the  results 
of  Kendrick  et  al.  (1977). 

The  results  of  the  binding  energy  calculations  indicate  also  that 
Li^so  should  be  unstable.  This  contradicts  the  experiaental  results 
of  Wu  (1976);  however,  other  recent  experiaental  works  indicate  that 
clusters  with  odd  nuaber  of  atoas  are  less  probable  than  those  with  even 
number  of  atoms.  The  origin  of  this  discrepancy  could  also  be  due  to  the 
paraaetrization  aethod  inherent  in  the  CNDO/BW  approximation,  as  correctly 
observed  by  Skala  (1981). 

The  Li^  results  are  probably  the  most  exhaustive  ab  initio  study  on 
this  system.  The  Cl  in  these  cases  were  fairly  large,  typically  11000 
configurations  counting  all  single  and  double  excitations.  At  the  Hartree- 
Fock-Configuration-Interaction  level,  the  binding  energy  increases  along 
the  following  sequence  of  geometries:  tetrahedral,  oblong,  triangular. 
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square,  linear  equidistant,  molecular,  T-shaped,  rhombic  and  parallelopiped. 
The  CNDO/BV  study  of  Sfcala  (1981)  with  the  same  parameters  found  the 
following  sequence  of  geometries:  tetrahedral,  triangular,  square, 
equidistant,  molecular,  linear,  oblong,  T-shaped,  rhombic  and 
parallelopiped.  Other  works  on  Li^  usually  take  into  account  only  two  or 
three  geometries  and  the  most  stable  system  is  that  of  square  fora  (Companion 
1978).  There  is  however  no  contradiction,  since  Li^r^°  or  Li^para 
were  not  considered.  The  binding  energy  2Etot  (L^)  -  Efcofc  (Li^para) 
is  found  to  be  0.495  eV,  to  be  compared  with  the  result  0.956  eV  of  Skala 
(1981)  and  0.668  eV  of  Companion  (1978). 

Finally,  the  result  of  the  dipole  moment  calculations  show  that  except 
for  Lij*so,  Li^mo^,  Li^  and  Li^para,  the  dipole  moment  equals  zero. 

Since  most  stable  Li^  systems  should  have  nonzero  dipole  moment,  experi¬ 
mental  determination  of  dipole  moments  can  help  determine  the  cluster 
geometries. 
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Table  1.  The  values  of  the  geometrical  parameters,  the  homo-lumo  gaps  AE, 
the  total  energy  EtQfc  in  a.u.,  the  binding  energy  per  atom  .in  a.u., 
and  the  dipole  moment  in  Debye. 


Cluster 

R 

AE 

E 

tot 

*b 

D 

U2 

2.672 

0.1903 

-14.7966 

0.0168 

0.0 

Lij1'1 

2.96 

0.0431 

-22.1654 

0.0070 

0.0 

Li3eq 

2.75 

0.0938 

-22.1696 

0.0084 

0.0 

Lij1'0 

R-2.82 

0-72.2° 

0.0674 

-22.1735 

0.0097 

0.375 

Li 

Lx4 

3.01 

0.1213 

-29.5773 

0.0128 

0.0 

U4tri 

3.89 

0.1218 

-29.5874 

0.0151 

0.0 

LiA8q 

2.84 

0.1223 

-29.5890 

0.0158 

0.0 

Li  eq 

Ll4 

2.65 

0.1609 

-29.5921 

0.0165 

0.0 

Li  ®Ql 

Li4 

2.73 

0.1601 

-29.5937 

0.0169 

0.037 

Li  lin 

Ll4 

R-2.734 

r-2.631 

0.1617 

-29.5897 

0.0159 

0.0 

3.23 

0.1221 

-29.5827 

0.0142 

0.0 

Li4T 

2.70 

0.1489 

-29.6021 

0.0190 

4.72 

Li  rho 

L14 

2.57 

0.1716 

-29.6048 

0.0197 

0.0 

Li  Para 
Li4 

Rl«R2"3.05 

0.1648 

-29.6114 

0.0214 

0.045 

a-50.4° 


*  » 
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Table  2.  Ionization  potentials  for  Lifl  clusters 


Cluster 

First  Ionization 

Potential 

Second  Ionization 

Potential 

Li2 

0.1799 

2.4258 

Lijtrl 

0.1237 

0.3633 

Li3eq 

0.1591 

0.3357 

Llj1*0 

0.1292 

0.3619 

u4‘*1 

0.112? 

0.2192 

u4tci 

0.1435 

0.2247 

Li/* 

0.1380 

0.2305 

u4*o 

0.1686 

0.2070 

Li™1 

0.1682 

0.2051 

u4lln 

0.1695 

0.2063 

U4obl 

0.1388 

0.2284 

LiJ 

0.1501 

0.2229 

Li  rho 

Li4 

0.1600 

0.2165 

0.1563 


0.2184 


